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Abstract 

This paper presents a complete algebraic analysis of the renormalizability of the d = 4 operator F^ v in the Gribov-Zwanziger 
(GZ) formalism as well as in the Refined Gribov-Zwanziger (RGZ) version. The GZ formalism offers a way to deal with gauge 
copies in the Landau gauge. We explicitly show that F^ mixes with other d = 4 gauge variant operators, and we determine the 
mixing matrix Z to all orders, thereby only using algebraic arguments. The mixing matrix allows us to uncover a renormalization 
group invariant including the operator F^ v . With this renormalization group invariant, we have paved the way for the study of 
the lightest scalar glueball in the GZ formalism. We discuss how the soft breaking of the BRST symmetry of the GZ action can 
influence the glueball correlation function. We expect non-trivial mass scales, inherent to the GZ approach, to enter the pole 
structure of this correlation function. 



1 Introduction 

QCD is the theory of strong interactions describing quarks and gluons which displays confinement at low energies. The 
mechanism behind confinement is still not successfully described. Even if one omits the quarks, the theory remains con- 
fining. Therefore, confinement is highly entangled with the dynamics of gluons, which makes glueballs very interesting 
objects to investigate. The existence of glueballs would be a pinnacle of the correctness of QCD, however, so far, there 
is still no clear experimental evidence for the existence of glueballs. This is mainly due to the mixing of glueball states 
with meson states which contain quarks. By increasing the statistics and/or by doing more involved experiments creating 
certain glueball states which cannot mix with quark states (oddballs), one hopes to uncover some clear evidence for glue- 
ball states. We mention a few experiments to demonstrate the general interest in glueballs: PANDA [1 1, BES III [2| and 
GlueX [3 1, ALICE at CERN @). 



The lack of experimental evidence has not stopped the community to widely investigate glueballs in various theoreti- 
cal models, see J5] and their references therein. Currently, theoretical estimates of e.g. masses of the different glueballs 
are compared to the lattice data. In lattice gauge theories, there is no doubt about the existence of glueballs and one 
can even work in pure Yang-Mills gauge theory [6 |. There are many phenomenological models which contribute to our 
intuition in glueballs. More direct contact with fundamental QCD can be made by identifying suitable gauge invariant 
operators, which carry the correct quantum numbers to create/annihilate particular glueball states Q- This is in accor- 
dance with the direct approach to study bound states in quantum field theory [8|. The mass of the glueball can then be 
determined by the leading singularity in its propagator which, if the glueball is stable, is just a simple pole. Of course, 
it is necessary to take into account non-perturbative effects, as glueballs are inherently connected to the non-perturbative 
region of QCD. One widely used method to estimate these propagators is based on QCD sumrules l|9j [TOj, while taking 
into account condensates, sometimes in combination with instanton or other nonperturbative effects. Also in holographic 
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descriptions of QCD, such glueball correlators have already been investigated, see for instance lfTTl[T2"l . 

In this paper, we shall concentrate on identifying a suitable composite operator , which is a renormalization group in- 
variant containing F^ v , representing the lightest scalar glueball. Let us explain how we shall take into account a particular 
source of non-perturbative effects. For this, we need a bit of background. As is well known, the Faddeev-Popov quantiza- 
tion of the Yang-Mills gauge theory was constructed in order to restrict the path integration only over gauge inequivalent 
fields. This restriction is translated at the level of the action by implementing a gauge, e.g. the Landau gauge d^A^ = 0, 
through the introduction of extra terms in the action, which in return break the local gauge invariance. In 1977, Gribov 
showed lfl3ll that this gauge fixing procedure in Yang-Mills gauge theories does not entirely restrict the path integration to 
gauge inequivalent fields, i.e. there are still multiple gauge copies which all fulfill the Landau gauge condition. More- 
over, it appeared that the infrared behavior of the gluon and the ghost propagator is strongly influenced when handling 
these copies. Therefore, there was a need for a formalism which took into account these Gribov copies, even if it would 
be only in a partial way. After a semiclassical treatment by Gribov in fPH , Zwanziger managed to construct an action 
which analytically implements the restriction to the Gribov region £2 fl4l . This action is called the Gribov-Zwanziger 
action Sgz- The region £2 is defined as the set of field configurations fulfilling the Landau gauge condition and for which 
the Faddeev-Popov operator, 

M" h = -d,(d,8" b + 8 r b A^ , (i) 

is strictly positive. Therefore, 

£2 = {A^d tl A a ^ = 0,M ab >0}. (2) 

The boundary, 3£2, of the region £2 is called the (first) Gribov horizon. The restriction of the path integral to £2 removes 
most of the Gribov copies in the Landau gauge related to (infinitesimal) gauge transformations [ 1 3 1 . However, there are 
still copies present in £2 and hence a further restriction to the Fundamental Modular Region (FMR), the region free of 
any Gribov copies, should be implemented. Unfortunately, till now, nobody knows how to handle such a restriction to the 
FMR. Therefore, the best analytical approach to restrict the number of gauge copies is by working with Sgz- We recall 
that Sgz is renormalizable to all orders lfl5l [161 [T7l . even in the presence of massless lfl8l [191 or massive quarks |20|. 
Implementing the restriction to the horizon introduces a first non-perturbative mass scale, the so-called Gribov parameter 
y 2 . Also, we have found in lF2Tll22 | that the auxiliary fields introduced by Zwanziger to construct the action Sgz, develop 
their own dynamics. This can introduce a second mass scale into the action. Generally, such non-perturbative mass scales 
are expected to be transmitted into the pole mass of the correlation functions. 

In a previous paper |23l we have investigated the operator Fi, in the ordinary Yang-Mills theory with Landau gauge 
fixing. This was already far from being trivial as at the quantum level mixing occurs with two other 4 dimensional oper- 
ators, i.e. a BRST exact operator £ = s(. . .), and an operator 9i which vanishes upon using the equations of motion. We 
have shown that this mixing does not have consequences when turning to physical states. Indeed, a BRST exact operator 
is always irrelevant at the level of physical states as the Yang-Mills action is invariant under the BRST symmetry. In 
this paper, we shall elaborate on the operator F^ v by investigating it in the more complex Gribov-Zwanziger framework, 
whereby exploiting the construction we have set up in [23 1. In this case, a similar mixing shall occur, but, in contrast with 
the Yang-Mills case this mixing shall have consequences at the physical level. Indeed, as the Gribov-Zwanziger action 
gives rise to a soft breaking of the BRST symmetry |22l . one can figure out that the corresponding BRST exact operator 
which will mix with J^,, will no longer be irrelevant. Let us mention that an attempt to calculate the glueball correlator 

•^vM^ap W) has been done in 11241 . but without taking into account the mixing of F^ v with other operators. We start the 
paper with an overview of the Gribov-Zwanziger action in section 2. We also recapitulate the Refined Gribov-Zwanziger 
action which takes into account the dynamics of the new fields introduced by Zwanziger. In section 3, a renormalizable 
action including the local, non-integrated operator F^ v (x) is constructed whereby in section 4 we shall analyze the mixing 
of this operator to all orders. In section 5, we shall determine the renormalization group invariant which contains F^ v . We 
end this paper with a conclusion in section 6, where we also present some insights on the potential relevance of the soft 
BRST symmetry breaking of the GZ action. 
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2 Overview of the (Refined) Gribov-Zwanziger action 
2.1 The original Gribov-Zwanziger action 

In this section we shall shortly recapitulate the ordinary Gribov-Zwanziger action in Euclidean space time which imple- 
ments the restriction of the path integral to the region £2. In [14], Zwanziger derived the following action, 

Sh = SYM+Sgf+Y 4 fd d xh(x), (3) 

with Sym the classical Yang-Mills action, 



Sym = l -Jd d xF« v F« v , (4) 



Sgf the Faddeev-Popov gauge fixing 



and h (x) the horizon function, 



The horizon condition: 



% = J d d x (b a d l A» + c?d fJ Dfc b ) , (5) 



h{x) = g 2 f bc A»(>M-Yf deC K- 



(h(x)) = d(N 2 -l), (7) 

with d the number of space-time dimensions, needs to be fulfilled in order to assure that we are working with a gauge 
theory quantized in the Landau gauge. This was proven using statistical arguments in lT4l[T5l . The action Sh contains a 

non-local term, but one can localize the horizon function by introducing the following set of additional fields: ^', c ,<p^ 

which is a pair of complex conjugate bosonic fields, and , co^ c ) , which is a pair of Grasmann fields. After this proce- 
dure, Sh gets replaced by Sgz, which reads 

Sgz = So+Sy, (8) 

with 

So — Sym + Sgf 

+ J x'xfadv (M a ;+g.f abn A b v y'; c ) -ca^v (9v<+g/^ m A*< c ) -g{d v m a ;)f" hm (D v c) b ^) , 

Sy = -fg J d d x (f abc Affi + f bc A% c + ^ (jy 2 - 1 ) y 2 ) , (9) 

We can further simplify the notation of the additional fields ^9^ e ,(p^ c ,ro^ c ,(B^'J as So displays a symmetry with respect 
to the composite index i = (fJ,c). Therefore, we can set 

(^,9^ C ><) = (W,<Pf,«M) , (10) 

so we get 

So = Sym + % + j d d x (cp?^ (z)fcpf) - c3f d„ (o>f ) - gf^d^fD^fj . (11) 



Finally, the horizon condition (0 can be written in a more practical version as 

dj 

whereby the quantum action F is obtained through the definition 



^ = 0, (12) 



e~ r = /[d4>]e- SGZ , (13) 
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where J [d<t>] stands for the integration over all the fields. 



For the Gribov-Zwanziger action, the conventional BRST symmetry is softly broken fl4l [22ll . We recall that the BRST 
transformations of all the fields are given by 

sA<; = -(D M cy , St? = l g ft>c c b c c j 

sc a = b a , sb" = 0, 

«Pf = <»>?, s<df = Q, 

sc5f = 9f, iff = 0. (14) 

The existence of this explicit breaking can be easily checked by releasing the BRST transformation s onto the action 

Sgz, 

sSgz = gfjd d xr bc (Ay;-(D^c m )(f fl c + i ? b ;)). (15) 
We refer to [22 1 for more details concerning this breaking. 

In order to discuss the renormalizability of Sgz, we treat the breaking as a composite operator to be introduced into the 
action by means of a suitable set of external sources. This procedure can be done in a BRST invariant way, by embedding 
Sqz into a larger action, namely 



£GZ = SyM + Sgf + So + Ss, (16) 



whereby 



d d x (-ufDftf - vfDfmf - u<;v f 
J d d x (-Mfrfftf - g f abc u^'D b l fc d v? c l + u£r>f co; 



~N M D p co, - V M D M <p,- +gf D M c co ; - -M^V p +U fJ N IJ j . (17) 
We have introduced 4 new sources t/° ! , V"\ and A^" with the following BRST transformations, and 

sU™=M™, sM^ = 0, 

sV*=Nf*, sN™=0. (18) 

This embedding into a larger action is necessary for the algebraic proof of the renormalizability as this heavily relies 
on having a BRST symmetry. Replacing the sources with their physical values in the end, returns the Gribov-Zwanziger 
action, 

= 7 2 8%v, (20) 

phys 



M ab 



yab 

phys 



as one can easily check. 



2.2 The Refined Gribov-Zwanziger action 

Let us explain the origin of the Refined Gribov-Zwanziger action. In the original Gribov-Zwanziger framework in 4 di- 
mensions, one obtains an infrared suppressed, positivity violating gluon propagator which tends towards zero for zero 
momentum and an infrared enhanced ghost propagator. This behavior of the gluon and the ghost propagator stemming 
from the action Sgz seemed to be in agreement with the lattice results for a long time. Until more recently, the authors of 
ll25l discovered a completely different behavior of the propagators in the deep infrared working on larger lattices. Now 
the ghost propagator no longer seems to be enhanced and the gluon propagator reaches a finite value at zero momentum. 
Since the publication of [25|, more lattice data have confirmed these striking results [26, 27, 28 29 30 31 1. Therefore, 
the Gribov-Zwanziger framework appeared to be in disagreement with these newest lattice data. However, in ||2T1|221 . we 
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have shown that it is still possible to obtain results with the help of the Gribov-Zwanziger action which are in qualitative 
concordance with these new lattice data by taking into account the dynamics of the fields ((p^ c ', (f^f, <3^ c ,<»° c ). This gives 
rise to additional non-perturbative effects within the Gribov-Zwanziger framework as, for instance, the dimension two 
condensate ((p^ e (p", c — S5^ C C0^ C ), which has been found fl2"Tll2"2l to be proportional to y 2 . It is apparent that the dynamics of 
these extra fields is highly entangled to the existence of the horizon. Therefore, we have refined the Gribov-Zwanziger ac- 
tion by explicitly adding the operator <p^' (f^f — (iojfco^ from the start, while preserving the renormalizability of the theory. 

The Refined Gribov-Zwanziger action is thus given by 

SrGZ = SGZ + S<p<p + Sen, (21) 

whereby 

S en = 2^fi- [d d xqfM 2 . (22) 

We have introduced a new parameter q and a new mass M 2 . The second term S en is a constant term, which is comparable 
with the term — y 2 J d d xd (N 2 — l) y 2 in the original Gribov-Zwanziger formulation (O. This term will allow us to remain 
inside the Gribov region £1. For more details on this construction, we refer the reader to l22ll . 



3 The (Refined) Gribov-Zwanziger action with the inclusion of the scalar glueball 
operator 

3.1 Generalities 

The most natural way to study the lightest scalar glueball is by determining the correlatoiQ ^ —j^- ^ • This correlator 

can be obtained by adding the operator F 2 V /4 to the (Refined) Gribov-Zwanziger action by coupling it to a source q(x). In 
this fashion, we obtain the correlator as follows, 

F\ X )F 2 { y)\ (23) 



8q(y) 5q(x) j q=0 

with Z c the generator of connected Green functions. In ll23l we have studied the glueball operator in the standard Yang- 
Mills theory, supplemented with the Landau gauge fixing. The framework we have set up for pure Yang-Mills theories, 
can be now extended to the more complex case of the Gribov-Zwanziger action, which is our current goal. 

Unfortunately, simply adding F 2 V to the action turns out to be too naive. In [23 1, we have seen that the 4 dimensional 
operator F 2 V mixes with other 4 dimensional operators in d = 4, in agreement with the general theory concerning the 
renormalization of gauge invariant operators l32l [33l l34l . 

Obviously, we also expect a similar mixing in the Gribov-Zwanziger framework. As outlined in ||23~1 [35l [36l . we can 
distinguish between 3 different classes of dimension 4 operators. The first class C\ is the set of the gauge invariant opera- 
tors, for example F 2 V . The cohomology of the nilpotent BRST symmetry generator s allows to identify the C\ operators 7 
as those which can be written as sf = 0, but also J ^ s(. . .). The second class C2 are the BRST exact operators, which 
are trivially BRST invariant due to the nilpotency of the BRST operator. Thus £ e C2 if and only if £ =*(...). The third 
class C3 contains operators which vanish when the equations of motion are invoked. One can then argue that the mixing 
matrix of these operators must be upper triangular, 

I Xrf yr Zy E Z^^ \ ( f \ 

(24) 

This particular behavior of the mixing of the various class of operators can be easily understood Il35ll36l . Bare C2 oper- 
ators cannot receive contributions from gauge invariant C\ operators: matrix elements of a bare BRST exact operator £ 





1. At least, this is our starting point. Later, we shall determine a renormalization group invariant %, containing F^, so we can calculate {5^. {x)%, (y)). 
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between physical states are zero. But, if there would be a renormalized gauge invariant C\ contribution in the expansion 
of £ , then there would be room for a nonvanishing contribution, which is of course a contradiction. Likewise, any C3 
operator vanishes upon using the equations of motion, while C\ - and a C2 operators in general do not, hence a C3 operator 
will not receive corrections from the other type of operators. 



In |f23ll , we have strictly proven in an algebraic fashion the upper triangular form of the mixing matrix for the operator F 2 y , 
just by using the Ward identities of the action. In particular, we have proven that the following action is renormalizable 
for ordinary Yang-Mills gauge theories in the Landau gauge, 



2-YMglue 



Sym ' 



d d xaA 



x (b'dpA" 

8(SYM + Sgf) 



—a *} T\ab b 
c o,.D,, c 



d d xkd ft c a A a p + J d d xr](d fl b a A^ + d ft ^'Dfc b 



ec 2 



5A« 



(25) 



ec 3 



whereby we see the three different classes of operators arising. We have introduced three new sources: the doublet (A,,r|) 
with sr\ = X and the color singlet a. The term (d^b"A^ + d p c"D' :d: 'c b ) is indeed an element belonging to the second 
class C2, as we can rewrite it as s{d fJ c a A" IJ ). In l23ll . we have introduced the last term through a shift of the gluon field 



A 1 ' 
A ,j 



oA£. 



3.2 Inclusion of the glueball operator in the Gribov-Zwanziger action 

With the mixing of the 4 dimensional operators in mind, we can propose an enlarged Gribov-Zwanziger action containing 
the glueball operator F 2 V . This action will turn out to be renormalizable. For this, we can make two observations. Firstly, the 
limit, {(p,(p, co, co, U,V,N,M} — > 0, has to lead to our original Yang-Mills action EYMglue with the addition of the glueball 
terms given by equation ([25l l. Secondly, setting all the terms related to the glueball term qF 2 equal to zero, we should 
recover the Gribov-Zwanziger action Eqz in equation ([Tol l. Therefore, we propose the following starting action: 



--glue 



£gz - 

+11 



d d xqF^ + J d d xs(r\ 



d d xqF» v F« v - 



d d x{X 



+ dwd<p + gf akh d(o a A k <p b 
d(dd(p + gf akh dm a A k q> b + 



v a D ab m b + uv 



U a D ab y b - 



-V a D" b T» b 



-UV 



M"K + d,c"Dfc b + acpacp - 3co3co + gf akb dT$,«A k q b + gf akb d® a D* a c a q>» - gf akb dWA^ 



\7K a A kft\b 



+ M M D p <P; +gU M f D^c (p f - Up Dp co,. +Np D M co,- -gV^ f D p c co ; + V M D M tp ( - 



abc r\ah b tl 



ab,J) 



tlb-rrb 



rabc r\bd d—c 



7 ai r\iib—b 



-M a 'V al -U m N m 



(26) 



Indeed, upon taking the limit {<p,(f>, CO, CO, U,V,N,M} — > 0, we recover the Yang-Mills actional (l25T> and setting all sources 
equal to zero (q, T|, X) — ► 0, we find our original Gribov-Zwanziger action back, see equation ( fl5b . Notice that in princi- 
ple, we could have taken other possible starting actions which also enjoy these two correct limits. We could have tried to 
couple different sources to the different BRST exact terms instead of employing only one source T|. However, this would 
not lead to a renormalizable action, while the action ( 1261 ) does turn out to be renormalizable, as we shall prove. 



We shall now try to establish the renormalizability of ( 126*1 ) by using the algebraic renormalization formalism |[37l . 

The first step is to introduce two auxiliary terms necessary for the process of renormalization. Firstly, we add an ad- 
ditional external term 5 ex t,i to the action, 



>ext,l 



J d d x (-K%Dfc b + ^gL a f abc c b A 



(27) 



which is needed to define the nonlinear BRST transformations of the gauge field A^ and of the ghost field c". K'J and L" 
are two new BRST invariant sources which shall be set equal to zero in the end, 



K\ t, =0, 

h Iphys 



^lphys = 0. 



(28) 



2. The term proportional to the equations of motion will be introduced later. 
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Therefore, these sources can be seen as two auxiliary sources which do not change the physics of the theory. Secondly, 
we also introduce the following external term, 

W - J &xs{&A!fi%) = Jd d x Y i A;dW;-Jd d x i K X I D; b c%W!+X I A a ^) , (29) 

whereby (X,-,!;) is a new doublet of sources, i.e. sX[ = Y,. This additional term is necessary in order to have a sufficient 
powerful set of Ward identities. Without this term, two Ward identities of the original Gribov-Zwanziger action would be 
broken which are absolutely indispensable for the proof a the renormalization of the action (see Ward identity 8. and 9. in 
the list below). Again, in the end, we shall set 

Xi\ phys = 0, Yi\ phys = 0, (30) 

We shall thus continue the analysis with the following action 

E = 2-glue + ^ext.l + ^ext,2 ■ (31) 



The second step is to search for all the Ward identities obeyed by the classical action E. Doing so, we find the following 
list of identities: 

1 . The Slavnov-Taylor idenitity : 

J(E)=0, (32) 

where 

f , ( S£ 8£ 8E 8E . 8E 8E . 8E • 8E ai 8E , 8E 8E \ 
5(E) = /^(s^SA^S^^^ (33) 

This identity is a functional translation of the BRST invariance s. 

2. The U(f) invariance: 

UijL = 0, (34) 

with 



-<af— -co" ^—r ~M"J— '—: -U, 



Using Qf = Uu, we can associate an extra quantum number to the /-valued fields and sources. One can find all 
quantum numbers in TABLEQ]and TABLE|2] 



3. The Landau gauge condition: 

8/7 



H = 9^-9^). (36) 



The modified antighost equation : 



5. The ghost Ward identity: 
with 



g a z = Ku (38) 





(39) 



7 



6. 



Two linearly broken local constraints: 



8E _8E_ 
Sip "SMtf 
SE 

"SM" 



SE 
Sco"' 



= 8f abc A b x i - ^\gf abc 4n i - WW® , 



(40) 



7. The exact %xj invariance: 
with 



^;E = 0, 



■89? 



" SM?< 8T> 



(41) 



An extra integrated Ward identity: 



5 of 5 



8X '8A, 5Z7" " 8M; 



(42) 



which expresses in functional form the BRST exactness of the operator coupled to X. 
The integrated Ward Identity: 

8 



J d d x I c a 



8co" 



?>K« ST, 



(43) 



10. The X-and T-Ward identities: 



/ 



d d x 



/, \ § ^§ — a § —a § 

(l-T 1 )_-X_ + C0 i ^ 7 + C P ,— 



d d x 



8Y 
(1-Tl) 



8c° 
8 



8 



, CO 

Si" ' 8b a 



0. 
0. 



(44) 





A" 
H 


c" 


c" 




<p? 




CO? 


CO? 


dimension 


1 





2 


2 


1 


i 


1 


1 


ghost number 





1 


-1 











1 


-1 


^/-charge 














1 


-1 


1 


-1 



Table 1: Quantum numbers of the fields. 











yen 


K 


L a 


9 


11 


I 


X 1 


Y' 


dimension 


2 


2 


2 


2 


3 


4 











1 


1 


ghost number 


-1 





1 





-1 


-2 








1 





1 


2/-charge 


-1 


-1 


1 


1 

















1 


1 



Table 2: Quantum numbers of the sources. 



Let us stress here that it is of paramount importance to have a good set of Ward identities to start from. For the construction 
of the action E, one should keep in mind the limits to the ordinary Gribov-Zwanziger case and to the Yang-Mills action 
with the inclusion of the glueball term. It is logical that an identity which plays a crucial role in one of the two limit 
cases, should not be broken by the action E, as E can be seen as an enlargement of the two limit cases. This is the reason 
why we have introduced S ex t,2- Without the auxiliary sources X{ and Yi, the extra integrated Ward identity d42l and the 
integrated Ward identity d43l are broken, and without these two identities one cannot prove the renormalizability of the 
action in an algebraic way. Let us also mention that in the ordinary Gribov-Zwanziger case, we have two extra linearly 
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broken constraints, belonging to the set of Ward identities in equation (f40b . However, it is not a problem that these two 
identities are broken, as the other two linearly broken constraints in equation d40b turn out to be equivalent at the level of 
the algebraic renormalization, namely: they have the same effect on the counterterm. 

Subsequently, we are ready to turn to quantum level. The third step is to characterize the most general integrated local 
counterterm E c which can be freely added to all orders of perturbation theory. E e is however restricted due to the exis- 
tence of the Ward identities. Let us investigate these restrictions a bit closer. The classical action changes under quantum 
corrections according to 



E^E + /zE c , 



(45) 



whereby h is the perturbation parameter. Demanding that the perturbed action (E + /zE e ) fulfills the same set of Ward 
identities obeyed by E, see ll37l . it follows that the counterterm E c is constrained by: 



The linearized Slavnov-Taylor identity: 

«rE c =0, 

where 2Jg is the nilpotent linearized Slavnov-Taylor operator, 

' 8E 8 8E 8 8E 8 SE 8 



d d xi 



dKfj 8A£ 8A£ dKjj 8L a be" 8c° bL a 



-MI 



" sua 



-n: 



-b a . 



,1 
ScSf 

,8 J' 



,1 
Sep? 



(46) 



and 



«£®£=0. 



(47) 



2. The U(f) invariance: 



(48) 



Utj is given in expression 



3. The Landau gauge condition 



8E C 



0. 



(49) 



4. The modified antighost equation: 



8E C 8E f / 8E f 



= 0. 



(50) 



5. The ghost Ward identity: 



with 



g a r = o, 



(51) 



(52) 



6. The linearly broken local constraints: 



8E' 
8^" 
8E £ 



- 8E e 
^8M? =0 ' 



(53) 



9 



7. 



The exact %ij symmetry: 



%E C = 0. 



(54) 



8. The extra integrated Ward identity: 




The integrated Ward Identity: 




10. The X-and F-Ward identities: 



/ 



d d x 



5 8 



CO 



/ 



d d x 



(i-n) 



5 




8 " 


8c fl " 




8 


h<3< 


8 " 


8i" " 


W 



E c 



(57) 



At this point, we are ready to determine the most general integrated local polynomial E c ' in the fields and external sources 
of dimension bounded by four and with zero ghost number, limited by the constraints d46l>— (T57T) . The linearized Slavnov- 
Taylor identity plays an important role in simplifying the form of the counterterm. Indeed, the counterterm can be param- 
eterized as follows: 

E e = («£ closed but not exact part) + SeA _1 , (58) 




whereby Ej is a cohomologically non-trivial part while Ej represents the cohomologically trivial part. A is the most 
general local polynomial with dimension 4 and ghost number — 1. One can prove that all fields and sources belonging to 
a doublet can only enter the cohomologically trivial part ll37l . This is exactly the reason why we have opted to introduce 
the source T], which is coupled to the BRST exact term, as part of a doublet. In this way, the source r| can only enter the 
trivial part, and turns out to be useful to explicitly prove the upper triangular form of the mixing matrix in equation (124-b . 
One can now check that the closed but not exact part is given by 

Ej = fl()SYM + V>YM, (59) 

whereby 

Sym = Jd d xq^F^, (60) 
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and the trivial part is given by the following rather lengthy expression: 



d d x 



ai (K"+ d^')A" + a 2 L"c" + a 3 U" d^f + a 4 V", 3„o3? + a 5 



+a,U;r; i + aigf bc U; i ^ + a %g f b X^ 

+a i iX'ffif 9A° + ai 2 X''ac^A° + cisX'ipfB" + a u gf abc T^p) + a\ 4 gf ahc X i <affi& c j 
+a l5 X , Wfb a + a l6 X i U , fJ "A^ + a ll gf abc X i lS? i tf<$ : j + a^gU.X'Wj^ + +a'l 7 gf abc X i ®yffi i 
+a x% gf ahc X i T3>c b c c + a 19 X i X i tfffi + a'^rx^W' + a^Y^M + a 2l gfabcY i &f& b j^ j 



+q 



b x {K« + d^)Al + cifd^ + b 2 L a c a + b 3 u;, drf + c^Uffl + b 4 V^d^ + c 4 d p Vffi 



+b 5 W'd 2 tf + csdjtildrf + d 5 d l O%tf + b 6 Uffi + b 7 gf bc U" p , q>?A£ + hgr%«a b A 



•abc j tci ,Jb a c 



rabC\TCl -rrrb \ C 



+b 9g f abc m^d^ b +c 9 gf bc ^{d,A^ b + 

+d 10 dX i W?A a M + b n X i <fic? + b l2 gf ahc X i (oftfm C j + b[ 2g f ahc X i (o i ;m b o? j + b l3 X i (%b a + buX'ufAl 

+b l5g f ahc x i my b ^ j + b^gf^TWj^j + b^gUcX^'qffi + b l6g f ahc x i ^c b c c +b 17 x i x i ^w j 
+b' 11 x i xj^m a j +b n x i Ym , M +b l9 gf ahc Y i w^ c j + b' l9 gf abc Y l w;w b ^ + b w Y 1 ^ 

+r\ eyK^Al + e'yd^'A^ + f^A", + e 2 L a c a + e 3 U% d^f + f&JJffi + e 4 V£ 9„c5? + f 4 d,V» CO? 
+e 5 C^d 2 cpf + ftd^Vi + gsd 2 ^^ + e 6 Ufift + e 7 gf bc U^ cp?A< + e^V^A^ 

+g xo dX i T5 c !Al + e n X'^c a + e 12 gf abc X i Oif CO^ + e' l2 gf ahc X i (d)m'm +e 13 X i O%b a + e 14 X 't^Aj 
+ei 5 gfabcX i O% <pi4 + e'^gfabcX^^ + e'( 5 gf abc X l Wf}^ + e l6 gf abc X i W!c b c c + enX'X^Wj 

+e' ll X i X%'Wj + ei & X i Y J W a W a + ei 9 gf a bcY i O%& -cpj + e' l9 gf ahc Y'm c jm b ^ + e^Y^T* 
igfabcXY^Wj + h' l gf abc X i <p ai ® l ffi + h 2 XV'la1 + h 3 mfm b (?^ b + (variants of h 3 ) 
The coefficients a,, aj, etc. are a priori free parameters. 



(61) 



As the attentive reader might have noticed, we did not include terms of the form (q 2 ...), (r| 2 ...), (qr\...), (q 3 ...), 
(kq 2 . . .) etc., into the counterterm. However, by just looking at the dimensionality, the ghost number and the constraints 
on the counterterm, one might conclude that certain terms of quadratic and higher order in the sources (q, T|, A,) are per- 
fectly allowed. One can imagine that an infinite tower of counterterms would then be generated and thence it would be 
impossible to prove the renormalizability of the action as new divergences are always being generated, which cannot be 
absorbed in terms already present in the classical action. However, we can give a simple argument why one may omit 
this class of terms with the help of an example. Assume that we would introduce the following term of order q 2 in the 
action, 



r F 2 



Subsequently, when calculating the correlator, this term would give rise to an extra contact term contribution, 

. 

term due to part in q term due to part in q 1 



5q(z) 5q(y) 



F2 tt F2 w\ +5(y _ z) /F 2 (y) 



(62) 



(63) 



Eventually, we are only interested in the correlator for z ^ y and therefore we can neglect the term (1621 ) quadratic in 
the source q. Moreover, when studying the case z = y, one should also couple a source to the novel composite operator 
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F 4 = F^ V F^, which is not in our current interest. We can repeat this argument for all the terms which are zero in the 
physical limit. Therefore, this argument is not only valid for the dimensionless sources q, r| and A,, but also for the massive 
sources K^, L p , Xj, Y,. Though, some care needs to be taken. Let us explain this again with an example. The modified 
antighost equation has the following form: 

5E C ' . 5E C / 5E C \ 

In this case, due to the term d^lp-, one compares terms of quadratic order in the sources ~ qK^ . . ., with terms of first 

order in the sources ~ q This identity can never be fulfilled is one immediately omits all terms of quadratic order in K". 

Therefore, we have chosen to keep all the possible combinations of higher order in the massive sources in the counterterm 
( |6H as there are only a finite number of combinations, while keeping in mind the higher order combinations of the di- 
mensionless sources. Only after imposing all the constraints, we can then safely neglect the terms quadratic in the sources. 

With the previous remark in mind, we can now impose all the constraints d48ll-(l57li on the counterterm, which is a very 
cumbersome job. We ultimately find 



- g f ahc v; a m%c c - gUcd^^'d^ 

+6, 1 d d xq L-^ + d^d^ + Kp^ + Mp^ Ufdrf+rfd^ + Vfdjtf 



+ay J d d xr\ (d^c" +M*dtf* - Ufd M (tf+NZ%e% + 

+ai J d d xX ^Jjp^ ai + VPfSO"' + ^ ^ ,03 ^ "'^ ^ ,9 ^ " + U^V^j -aij d d x ^X'd^d^iA . (65) 

Only now, we can discard the term ~ qK^df,c a as it is of quadratic order in the sources. One could argue that we can 
also neglect terms of higher order in U™ and N° l . However, both sources belong to a BRST doublet. Moreover, the 
corresponding partner sources, M^ l ,V^ 1 , acquire a nonzero value in the physical limit, and it would be impossible to write 
the BRST exact term in our starting action E g i ue (see expression d26i i) as an s-variation when neglecting these kind of 
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,85 



,55 



YM 



terms. In summary, the expression 

E c = aoSYU + hSYM + ar I d d x +d l ^d fi c" +K a ^c" + M«'d^ - ^%G>™ 

+dj?%tf + ax%< + v*m* - ufrf - 8 fabcUfa%c c - gf ab cv; a m'%c" - gu^^%c c 



+ ay I d d xX^u;% l ^Uvp,w i, +d^%^ , +u; i v; i j - ai J d d * 

gives the general counterterm compatible with all Ward identities. 



(66) 



We still need to introduce the operators belonging to the class C3, which are related to the equations of motion, see 
section [XT] Therefore, the next step is to perform a linear shift on the gluon field A"^ in the action E 

A^A^ + oA;, (67) 

whereby a is a dimensionless new source. As this shift corresponds to a redefinition of the gluon field it has to be 
consistently done in the starting action as well as in the counterterm. Later on, we shall see that introducing the relevant 
gluon equation of motion operator through this shift, will allow us to uncover the finiteness of this kind of operator. 
Performing the shift in the classical action yields the following shifted action E' 



E' = S YM + / (b a d^ + Z%Df c*) + J d d x (-1% {D^cf + l -gL a f" K c b c 
+ j d d x (9? 3yDf cpf - cof 3 v £>f oof - gd^f bm D b v d c d ^f) 
+ J d d x{-Mfiy^^ b i -gU^ i f^ c D bd cS c i +U^ i Lif(Si b i 



-NfDfm\ - v; i D c ' h tf+gvff ahc D b ; l c d m < ; -m*v* + u a ^ 



d d xqF° v F° v 



d d xX 



d^A" + 3©3cp + gf akb dWA k (f> b + U"D ab y b + V a D ab o5 b + UV 



d d xr\ 



M"K + d^D" h c h + 3<p3cp - 3©3co + gfakbdfAV + gf akb dm"D kd c d y b - gf akh dWA k a 



■a Ak fn b 



-rClrJid <d r Jb 



K a A k M b 



+M"J (D^tY + gU"T bc (D M c) b tf - U a J (D^T +N a J (D,fik) a - gVff bc {D,c) b W, + V" 1 (D^,)" 



-M C "V C " - U"'N C " 



d a x [YiA a ^ -XiDfc'd^! +X i A«d ll <tf 



d d xaA" 



55 



YM 



d d xa 



d d xa 



d d xa 



d d xa { -d,b"Al + gf akh A k /d^ } 
-gfattdjftAtf +gf akb d fl G%Ay - g 2 f abm f b kA& a <? m Ay 



g 2 f ab cf b kdUWA k c d + g 2 f abc f bkd V a m c A k c d 



(68) 
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Notice that we have neglected again higher order terms in the sources ~ (ocr| . . .), ~ (aX . . .) and ~ (aq . . .) as the argument 
d63l is still valid. The corresponding counterterm reads: 



+b, I d d xq +d^d,c a +M«%rf u*d„a# +a^< + v;%<?; 
+ai J d d xr\ (d^c" +M a ;d^ - u;%<tf +<a A ,< + v£%i?j 

+^(j>%< + + V*M* - « - gfabcUfa%C C - gf abc V^'W%C C - gfabc^^C 

+ ai J d d xx(ufd^ ai + V^%^' + d^%(? ai + U^V^ -aij d d x(x%o?%c a j 



-ao / d x\ aA - 



+a , J d d xa (lA^d^ - 2A ( 'd 2 A; + 9gf abc A<ffid^ + 4g 2 f abc f cde A"A b v 4A^j , (69) 
once more dropping higher order terms in the sources. 

The final step in the renormalization procedure is to reabsorb the counterterm E' £ ' into the original action E', 

Z(g,<o,hQ) +hYF = EUo,«>o,*o,*o) + 0(h 2 ) , (70) 
We set = (A°, c", c* b a , tpf , CO?, f?, C5f) and 4> = (K^, L", M°', N™, Vf, Uf, X) and we define 

go = Z g g, 0O=ZJ /2< I ) . *o-Z4,4>, (71) 

while for the other sources we propose the following mixing matrix 

qo \ ( Z qq Zqr\ Z qJ \ / q 

110 = \ Z n y T] | . (72) 

■^0 / \ Z Jq Zjt| Zjj J \ J 

If we try to absorb the counterterm into the original action, we easily find, 

z g = i-h ao 



zf = l+h 



2 ' 



(y+fli), (73) 



and 



zi/ 2 = Z y 2 =z-^z^ 2 =i-h^ 



Zb — z A 1 , 
Zk — Z l J 2 , 

Zl = Zf, (74) 
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The results ( |73T l are already known from the renormalization of the original Yang-Mills action in the Landau gauge. 
Further, we also obtain 

7 l/2 1/2 -1/2 -l/4_ 1 , ai 

— ^ ~ .? A A — 1 ""2" ' 

1/2 _ -1/2 
7I/2 _ 7 -l 

Z M - 1 - — - Z g z A 

7 _ 7-I/2 

v 1 ; fll -7-1/2^-1/4 

Z v = 1 -h— = Z g ' Z A ' , (75) 

which are known from the original Gribov-Zwanziger action, see [15|. In addition, we also find the following mixing 
matrix 

l+h(b -ao) 

hbi 1 I , (76) 

hbi 1 



Zqq 


ZqT) 


Z q j 


Zx\ q 


Zr\T] 


Zy\j 


Zjq 


Zjr\ 


Zjj 



-1/2 -1/2 1/2-1/4 

Z\ = Z c Z A = Z g Z A . (77) 



while for the Z-factor of X we have 



Also this part was already known, see 11231 . So far, we have proven that the two limit cases are at least correct. Finally, we 
find the new results 

Zy = Z g Z A 12 , 

Z X = Z l J 2 Z A l ' A . (78) 

In summary, the action Y! is renormalizable. Moreover, we have only 4 arbitrary parameters, ao, au bo, b\, which is the 
same number as in the limit case {cp,cp, CO, CO, U,V,N,M} — > 0, i.e. the Yang-Mills case with the introduction of the glueball 
operator ~ F^ v ||231 . This is already a remarkable fact. 

3.3 Inclusion of the glueball operator in the Refined Gribov-Zwanziger action 

In analogy with ||22][2T| we shall add the two dimensional mass term ~ ((p?cp" — cofco-') to the action E g i ue in equation 



2-Rglue — 2Tlgi ue + + E en , (79) 

whereby 

L m = J d d x (s(-mW)) = I d"x(-J (cpfcp? - C5fa>f)) , 

£ en = [d d xqej, (80) 



with / and 9 new sources, and q the parameter already defined in equation ( f22l ). In order to agree with the physical action 
(EH . we define the following physical limit, 



We further define sJ = and s® = 0, hence the BRST invariance is guaranteed. 
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Let us now investigate the renormalizability of action ER 2 i ue . We can go through the same steps as in the previous sec- 
tion. Therefore, we again add the two external pieces, 5 ex t,i and S ext ,2 as defined in equation d27b and d29l l, to the action 

ERglue 

= ^-Rglue + ^ext, 1 + ^ext,2 ■ 

(82) 

Subsequently, one can easily check that all Ward identities ( f34b - ( f4TT > and (l44t remain unchanged up to potential harmless 
linear breaking terms. Therefore, the constraints d48l - d54b and d57b remain valid. Unfortunately, the extra integrated 
Ward identity d42l) and the integrated Ward identity (03]) are broken due to the introduction of the mass term. However, the 
mass term we have added is not a new interaction as it is only quadratic in the fields. Therefore, it cannot introduce new 
divergences to the massless theory E, and it can only influence its own renormalizatiorQ as well as potentially vacuum 
terms, i.e. pure source terms. Also, next to Ward identities ( f34b - (|4TT > and d44b . we have a new identity 

H - 

which is translated to the following constraint at the level of the counterterm, 

§ - 0. (84, 

As a consequence, L C R is independent from the source 0. Therefore, it follows that the form of the counterterm E^ can be 
written as 



E^ = E c + E}, (85) 

whereby E c is the counterterm (l66l of E and T,j is depending on /. One can now easily check that Ifj = kJ 2 , with K a new 
parameter as this is the only possible combination with the source J, which does not break the constraints d46b - d54l l and 
d57j. 

K is in fact a redundant parameter, as no divergences in J 2 will occur, as explained in 11221 . Therefore, the counterterm E^ 
is actually equal to E c . Defining 

Jo=ZjJ, (86) 

we find 

Zj=Z-' =Z g Z l /\ (87) 
and we have proven the renormalizability of the action E c/ . 



4 The operator mixing matrix to all orders 
4.1 Preliminaries 

Let us return to the mixing matrix of the sources q, r| and J and pass to the the corresponding operators. We have found 
that 

/ q \ ( Z qq \ / q \ 

lo = \ Z Jq 1 T( . (88) 
r J \Z Jq 1 J \J J 

We shall further need the inverse of this matrix, 









f z 


i)- 




Zqq 




Z J'I 
















j qo 


1 





no 





1 J 


V Jo 



(89) 



3. We employ massless renormalization schemes. 
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We can write the final action E' from equation d68l l in a more condensed form as 

£' = E G Z + Sext. 1 + 5 ext,2 + J d d 'x (qf + T| £ + OUf{) + J d" ' xkVt , 

whereby we have defined the operators 



T 

•E 



A A- Aa ' 



(90) 



(91) 



with 



9C 



d M c"A 1 ' + 3co8(p + gf akb dW'A k y b + U a D ab <p b + V"D ab o5 b + UV 



ar\ab w Jb i t txx r\cibyrh 



(92) 



It is then an easy task to construct the corresponding mixing matrix for the operators themselves. We recall that insertions 
of an operator can be obtained by taking derivatives of the generating functional Z c (q,T[,J) w.r.t. to the appropriate source. 
For example, 



To 



8Z c (q,%J) Sq 5Z c (q,r\,J) 5r) 8Z c (q,r\,J) SJ 8Z c (q,r\,J) 



bqa 



Sqo Sr) 



Sqo 



SJ 



and thus 



To 



Zqq 



(93) 



(94) 



and similarly for £q and M®. Henceforth, we find 




z- 1 
m 






-1 



ZjqZqq 
1 






(95) 



This is a nice result as we recover the expected upper triangular form. In addition, as £ has a Z-factor equal to 1, we 
also find that the BRST exact operator £ does not mix with H , although this mixing would in principle be allowed. This 
can be understood as follows. The integrated BRST exact operator £ is in fact proportional to a sum of four (integrated) 
equations of motion terms and two other terms, 



J d 4 x [d^b a A; + d,c a Dfc b + 3cp8cp - 3co3co + gfaufflA'tf + gf akb dm a D kd c\ b - gf akb da?A k <a b + irfrftf 



j jai rabc r\db Jb rf JC 7 jai r\ttb M h 1 \rai r\(ib—b T T ai rabc r\bd d-rrc , t rai r\ab—b , i/OIt T ai j jai\jai 

-g u n f D M c <P; - U M D M co ( - + N M D M co, - gV M f D M c co, + V M D M cp,. +M /J V^ I - N M 



d 4 x b a 



SE, 



GZ 



,8E 



GZ 



Sb a 



8c" 



— a §^GZ , —a §^GZ , ,, a i^GZ , rrC „§^GZ 



sm a 



" 8M" 



" sua 



(96) 



and therefore, like H , it does not mix with the other operators. Notice that we can rewrite the integrated BRST operator 
in two other forms: 



(HI 



■/*»(. 



Sb« 



8c" 



Sep" 



8co" 



" SV" 



(97) 



.flSEcz . a 8^cZ _„8Egz . „§£gz . .,„,§£gz . Art „ 5^Gz\ 
b — + c — +cp +co - +Af„ , ; +AC 1 



Sb" 8c" 



8co" 



5A^ J 



(98) 



Remark 

We can also use the refined action Ergz instead of Egz- We define Ergz as 



Ergz = Egz + Emm + E, 



(99) 
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whereby E^> and E en are defined in equation (fSOb . Replacing Eqz by Ergz does not alter equation (|95k but it does slightly 
modify expression d%l l, 

/" ,4 /"j4 / , a S^RGZ , _^§Ergz , ^,oE R gz . —a §Ergz . ..„,8Er G z , 7rfl ,5E RG z §Ergz , „§Ergz\ 

(100) 

and analogously for expression d97l i and ( |98l l. 



4.2 The physical limit 

In the next subsection, we shall work in the physical limit as our final intention is to examine n-point functions with the 
(Refined) Gribov-Zwanziger action itself. In the physical limit, £ becomes: 

,t* I ^ jajfl i ^ —ar\ab b , 21 7^ar\ub, n b ~\ -rrar^ab.^b , jrabc~\ -rrcir^bd d ,_c i .2 rabc \a.Jjc , .2 rabc Ad-^rbc 

£ l P hys = <V> A^ + dpC c + <Vp,-D /I (p ; . -d^co, ^ CD, + g/ d^D^ c cp ( - +r^/ A^ +Ygf A^ 

+ d(N 2 -l)y 4 . (101) 

From this point, we can omit the constant term d (N 2 — l) T 4 as it shall not play a role in the calculation of the glueball 
correlator. Later, we shall determine the renormalization group invariant %.{x) which contains F 2 v (x). As £ mixes with 
F 2 v (x), this renormalization group invariant shall also contain this constant term. However, a constant term can never 
contribute to the final glueball correlator (x)%_ (y)) as it can never help to produce connected diagrams between the two 
space time points x and y. Therefore, we shall simplify the calculations by omitting this term already from this point. 



In the physical limit y{ is given by 



whereby Sgz is the physical Gribov-Zwanziger action ([8]). Naturally, the mixing matrix d95l l stays valid. 



4.3 The mixing matrix to all orders 

It this section, we shall determine the mixing matrix d95l l to all orders. This proof is very elegant as it does not require 
to calculate any loop diagrams, and it is purely based on algebraic manipulations. We shall extend the proof given in 
1 23 1, which is based on ll38l . Moreover, as a byproduct, the proof shall also reveal some identities between the anomalous 
dimensions of the different fields, which can serve as a check on relations as in d74l ) and d75l) . We shall directly work with 
the physical action Sgz- in the end, we shall also look at the Refined Gribov-Zwanziger action, Srgz- 

We start our analysis with the following generic n-points function 

g"(x u ...,x n ) = (to(xi)...ty(xn)) = f[tyMxi)...ty(xn)er Sa * 1 (103) 

whereby fc, z = 1...8 stands for one of the eight fields (A» c a , c a , b", <pf, (i)f, <pf, Tflf), i.e <|>i = A,,, . . ., <|> 8 = Tnf . We 
shall immediately omit the vacuum term "f{N 2 — l)d in the action Sgz, as it is relevant only for the calculation of the 
vacuum energy and not for the calculation of n-points functions. The total number of fields is given by n, 

8 

n = J>, (104) 

i 

with rii the number of fields (]), present in the n-points function ( 11031 >. We are therefore considering the path integral for 
a random combination of fields. Subsequently, from the definition ( 11031 >, we can immediately write down the connection 
between the renormalized Green function and the bare Green function, which is, in a very condensed notation, 

S n = Il\ ni/2 So- (105) 

1=1 

From the previous equation, we shall be able to fix all the matrix elements of expression (|95l l, based on the knowledge 

AG" 

that -j-y must be finite in a renormalized theory. 
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We shall therefore calculate this quantity. The first step is to apply the chain rule: 



8 / dZ, 

E 

7=1 



-nj/2 



Hj/2 



d 8 2 dgl d 8 2 df a 



So" 



Next, we need to calculate the derivatives w.r.t. g 2 . 

• Firstly, we need to find dg^/dg 2 . We employ dimensional regularization, with d — 4 — e. If we derive 



«3 



£ 7 2 2 

A* z„g , 



w.r.t. /j and g 2 , combine these two equations and employ the following definition of the p-funtior0 

dg 2 



we obtain 



dg 2 



-eg 2 + p(s 2 



Secondly, we calculate We start from 



To = Z^y 2 



whereby Z^i = Zy = Zm due to the limit ( fT9l ). Deriving this equation w.r.t. g 2 yields 
dy 2 dZ 2 2 d\nZ 2 , 1 3^ 3lnZ y 



dg 2 dg 2 r ~ dg 2 To ~ n dg 2 ** dfi r ° 



and we have defined the anomalous dimension of y as 



§ y 2 



3lnZ 



-zg 2 mg 2 ) 



VYo> 



2i y 



Finally, we search for "^ 2 /dg 2 . Applying the chain rule gives 



dZ 



-n/2 



dg 2 



n 



z:r 12 dzl! 2 



1/2 3^2 



n ^z l J 2 

-Pi/2 4>' 
3^2 



fllnZ 



1/2 



Next, we derive — from the definition of the anomalous dimension, 



3lnZ 



1/2 



Yd,' 



3/j 3/j dg 2 

From expression dl 131 ) and ( II 14b . it now follows 

az: p ' /2 



a„2 8lnZ' /2 5lnZ'/ 2 

* * = Hs 2 +P(s 2 ))- * 



3g2 



3g2 



= -Pi 2 ,, 



n /2 y 



-zg 2 + Mg 2 ) 



4. We have immediately extracted the part in £. 



Inserting equation ( 1 109b and (11 15b into expression ( 1106b , we find: 

,-«i/2 / „ 



dg" 
d^ 



n,-z, 



(116) 



The right hand side still contains bare and therefore divergent quantities. We would like to rewrite all these quantities in 
terms of finite quantities so that we can use the finiteness of the left hand side to make observations on the right hand side. 
Also, we should rewrite in some manner the number n,- as the mixing matrix ( 195b is obviously independent from these 
arbitrary numbers. 



Therefore, as a second step, we shall rewrite the right hand side of (II 16b in terms of a renormalized quantity. Firstly, 

we calculate -X^Qn - Using 

3go 



1 (Ffo) 



1 (. . . Sfez 



«oCv) 



-cooCy) 



§5, 



GZ 



8coo(y) 



(117) 



we can write, 



2^0 

So~, — 



i 



1 I „ { — 5.V, 



GZ 



ScooM J 

§^GZ 



We have introduced a shorthand notation for an insertion in the «-points function, e.g. 



n (F 2 (y) 



V(xi)...V(z„) 



4 | \ 4 

Secondly, we analogously find 

_a_ 

Thirdly, we rewrite njQ^ by inserting the corresponding counting operatofl into the Green function, 



fi^gS = /d 4 y(^ "{^or & ^,o< + ^or &c A^o} 



§Sgz 



6^0) 



(118) 



(119) 



(120) 



(121) 



Inserting (II 18b . ( 1120b and (1121b into our main expression ( 1116b results in 



dg" 
~dg T 



-es 2 +l% 2 ). 



d J y 



j=i I 5%M ■ 



s" «>o(y) 



_55gz_1 e 



§Sgz 1 



'5^oWj 2 y r uv; 8coo(y)/ 2 9 t~ uw 5co M J 



(122) 



Notice that we have also absorbed the factor Yli^, 2 i nto me Green functions, and therefore we can replace gfi again 
by Q n . Finally, we need to rewrite all the inserted operators in the zi -points function g" in terms of their renormalized 
counterparts. For this we return to the mixing matrix d95"b and parameterize it as follows 



(123) 





5. It is easily checked that f d*y$ J -~ counts the number of <|)q insertions. 
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Here we have displayed the fact that the entries associated with a(g 2 ,e) and b(g 2 ,e), which represent a formal power 
series in g 2 , must at least have a simple pole in 8. Therefore, we can rewrite 

= (-e- a )T(y)+b<E(y)\ phys + bA(y)? Sl:Z 



-ejFoCy) 

I phys 



5AM 



My) 



My) 



8Sgz 
8A(y) ' 



(124) 



whereby we recall that we are working in the physical limit and we have replaced H | phys by the expression ( 11021 i. 
Subsequently, 



(125) 



as one can check with the Z-factors in d75l l. Finally, all the other operators are equations of motion terms, which appear in 
expression ( |96] l, d97| > and d98l and therefore have the same Z-factor as the operator <E , i.e. Z = 1 . Summarizing, expression 
(11221> becomes: 
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+ V £ " {yV/^A^ + fgr^Affi } 
where we have immediately taken the full expression of £ | p h ys in equation ( llOU . 



(126) 



From expression (11261 ). we can determine a(g 2 ,e) and b{g 2 ,s). As ^ T is a finite expression, we know that the right 



hand side of equation ( 1126t must also be finite. Therefore, as all the Green functions are expressed in terms of finite 
quantities, we can choose a set of linearly independent terms and demand that their coefficients are finite: 

5S G zl . e/2 + b-y A {g 2 ) 
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-Yip -8 y2 -b 
~^g T +W)' 



We can rewrite the coefficients of g" {<J } and g" {A^jjp j in (1127al > as 



—e — a 



1 (1+a/e) 



e/2 + /?-y A (s 2 ) 1 l+2(£-y A (g 2 ))/e 



-eg 2 + P(s 2 ) s 2 l~P(s 2 )/(eg 2 )' 
Hence, in order to be finite, we must conclude that 



-e# 2 + p(g 2 



2g 2 l-V(g 2 )/( £ g 2 



(127a) 
(127b) 
(127c) 
(127d) 
(127e) 

(128) 



a(g 2 ,e) 



P(s 2 ) 



,e) = TSi(r)-o 



1 P(s 2 ) 



2 S 2 



(129) 
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Notice that a and b depends on g 2 , but not on e. Therefore, the matrix elements of the first row of the parametrization 
(|140t only display a simple pole in e. 



Moreover, from the other equations we shall obtain relations between the anomalous dimensions of the fields and sources. 
Let us start with the coefficient of Q n {bd^A^} in equation ( |127bt , yielding 

-e/2-b-y b (g 2 ) 1 l + 2(b + y h (g 2 ))/e 



-«&2 + P<g*) 2g 2 l-V(g 2 )/(£g 2 ) 

which means that 



lp(^ 
2 g 2 



2 ,e) = -Y^ 2 )-^^- (13D 



Inserting the value of b{g 2 1 €) from expression ( 11291 > gives the following relation 

lA+lb = 0. (132) 

This relation is a translation of the relation Z l J 2 Z l J 2 — 1 found in equation (f74"l >. Indeed, deriving both sides w.r.t. /j 
gives 

-J-^A (zf zf) =JA+J b = 0. (133) 

Z A L b 

Analogously, for the coefficient of Q n {c^Z^fc*}, we find 

b{g 2 ,£) = -Yc-Yc, (134) 



yielding 



Ya + Yc+Yc = A, (135) 
2g z 



which is a translation of Z^ 2 ^ 2 Z^ 2 Zg = 1 as fi^L = — ^j. Next, the coefficients of d!27cb and d!27db lead to 

P P P 

Y<p + Y<p+Ya = ^2 , Ja + Jm + jA = ' Yc + Yw + Yv + Ya = pr, (136) 

stemming from 

1/2 1/2 1/2 1 1/2 1/2 1/2 . 1/2 1/2 1/2 1/2 . ™ 

^9 A A A £ — 1 ' ^co A A A?" 1 ' ^ ^<P ^A A? — 1 ■ l."/; 

These relations originate from the relations derived in d74l i and (l75l l. Finally, the coefficients in equation d!27et are finite 
if 

- ■ ft..- = - V,„ - (Y:. : h Vi I y" i ^ 



Y<p-8 y 2 = -Ycp-5 f = ft = y 4 (/)-I^1, (138) 



or equivalently 



zV 2 z)[ 2 Zfa = 1 , z; /2 zj /2 Z,Z Y 2 - 1 , (139) 



which is also fulfilled as Z^ = Zy = Z g 1 2 Z A '^ 4 . 

In summary, we have determined to all orders the mixing matrix d95l l. For notational simplicity, we take the value (11 34b 
for b and we use the equality y c = y-: 



p(i' 2 ) 


2Y, 




Eg 2 
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1 



(140) 

We have encountered numerous checks which show the consistency of our results. 
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Remark 

This matrix is also valid for the refined action 5rgz- One can repeat the proof by replacing 5gz with 5rgz and by adding 
the following term in M 2 — J to the game, 



-M 2 J d d Jc(^fcpf — cSfcof) , 
see equation ( 1221 ). In the end, expression ( 1126b will collect an extra term 
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where we have introduced the anomalous dimension of M 2 , 

§ M2 = 



3lnZ 



M- 



This leads to the following extra coefficients 



g"{-M 2 ^} 



so that 



or equivalently 



y«p + y<p + 5 M 2 =o, 



zi /2 z l/2 z , - 1 

_ 7 7 V2 



£"{M 2 co>?} 



-£S 2 +P(s 2 ) 



Yoj + Yq) + 5 M 2 =0, 



1/2 1/2 1/2 
(0 A B> M 2 



which is correct as Zj = Z M i = Z g Z^ , see equation (|87| |. All the other relations stay valid of course. 



860) 
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(144) 



(145) 



(146) 



5 The glueball correlator 

5.1 A renormalization group invariant 

As the final step of our analysis, we shall try to determine a renormalization group invariant operator which contains 

f =^p. This is useful as we would want to obtain a renormalization group invariant estimate for the the glueball mass, 
i.e. the pole of the corresponding correlator. This analysis is completely similar to the one presented in l23l . due to the 
fact that the mixing matrix Z is exactly the same. However, for the benefit of the reader, let us repeat the analysis. We 
define the anomalous dimension matrix Y of the mixing matrix Z as 



(147) 



With the following derivatives, 
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(148) 
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we obtain 



/ „2 3(P/g 2 ) _7„2 3y c _7„2 3T £ . 





(149) 



V 



o 







Notice that this anomalous dimension matrix is finite, as it should be. This matrix T is related to the anomalous dimension 
of the operators, since 

X = ZX => = /j^-X+Z/j^- 

Ofi Ofl 



dX 

=> (ij- = -TX, 



with 



X 



'J 

■E 

H 



X 



To 



We now have all the ingredients at our disposal to determine a renormalization group invariant operator. We set 



with k, I and m functions of g 2 , to be chosen in such a way that 



a dk , 2 d($/ g 2 ) 

fij-H = fi—T -kg 2 ^| 
dfi afi og L 



, dy c d£ dm 

og A ofj afi 



T + 2kg L ^<E +2kg A -^X +fi— <E +M— X = 0, 



hence 



£ = m . 



(150) 



(151) 



(152) 



(153) 



We therefore choose 



and we conclude that 



% 2 ) = ^, 

J(g 2 )=m(g 2 ) = -2y c (g 2 ) : 

Ms 2 ) 



-7 -2y e (g 2 )£-2y t (gV 



(154) 



is a renormalization group invariant scalar operator containing F 2 V , in the case of the Gribov-Zwanziger action Egz as 
well as in the case of the refined action Ergz- 



5.2 Irrelevance of the terms proportional to the equations of motion 

As we have found a renormalization group invariant, the final goal [39] shall be that of evaluating the glueball correla- 
tor 

GOW = (^(^p-f(x)-2y c (g 2 )'E(x)-2y c (g 2 )! f {(x)j (^fW-2y c k 2 )£W-2y [ .(gV(y))|> , 

(155) 

using the (Refined) Gribov-Zwanziger action. However, this is beyond the scope of the present article as this calculation 
shall be far from trivial, even at lowest order. 
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As usual the equation of motion terms like !tf will not play a role. Let us demonstrate this with a simple example, 

= ...8(jt-y) + 8(0)(j(*)), (156) 



8A«(y) ( 



which is zero as x 7^ y and 5(0) = in dimensional regularization. Therefore, expression d 155b reduces to, 

<^(*H(y)) phys = H&j (jW^(y)) + (2y e (^)) 2 ( £ (x)£(y)) phys 

-2y t (g 2 )^((j(x)£(y)) phys + (£(x)j(3;)) phys ) . (157) 

6 Summary and discussion of the relevance of the soft BRST breaking 

F 2 

In this paper, we have scrutinized the glueball operator J = -j- using the (Refined) Gribov-Zwanziger action Sqz (Srgz)- 
For this, we have followed the framework of an earlier work ll23l where we have investigated this operator for the more 
simple case of the usual Yang-Mills gauge theory, quantized in the Landau gauge. However, this framework is heavily 
based on the existence of the BRST symmetry while neither Sgz nor Srgz are BRST invariant ll22ll . Therefore, through- 
out the paper, we have relied on the extended model Egz and Ergz- With these "enlarged" actions, one can then draw 
very similar conclusions as in the ordinary Yang-Mills case. The results of interest, i.e. those for the (Refined) Gribov- 
Zwanziger action, then easily follow from these extended models in the physical limit, in which case certain external 
sources are assigned a suitable value. 

Firstly, the classically gauge invariant operator F 2 V mixes with two other operators, a BRST exact operator, £ = s [d^'A^ + 
3co3cp + gf a i ( bd&"A k ^ b + U"D ab y b + V"D ab & b + UV], and an operator proportional to the gluon equation of motion, 

using the algebraic renormalization procedure, we have determined the form of the mixing 

matrix Z to all orders, 

Zqq ~Zj q Z q q —ZjqZ q q 

1 I I L I , (158) 

1 





which has an upper triangular form, as required 

In a second part of the paper, we have completely fixed all the elements of this mixing matrix, by using only algebraic 
arguments. We have found 

, _ P(.? 2 ) 2y c (g 2 ) 2y c (g 2 
Eg 2 e 6 

1 ! - ; i- S9 > 
1 

which is completely analogous as in the case of the ordinary Yang-Mills theory ||23ll . This is already a remarkable fact. 
In addition, we have also encountered numerous checks on our results as we have recovered multiple known relations 
between the anomalous dimensions of all the fields and sources. 

In the final part, we have determined a renormalization group invariant including F^ v , given by 



m = ^ r L T-2y c (g 2 )<E-2y c (g 2 )?{, (160) 
8 



which is the main result of this paper. This operator would then be a good point to start the study of the (lightest) scalar 
glueball from, by means of the correlator (li (x)^K.(y)) pilys (39). 
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In standard Yang-Mills gauge theories, gauge invariant operators J only mix with BRST exact and equation of mo- 
tion type terms. While the latter always yield trivial information at the level of correlators, the BRST exact pieces drop 
out due to the BRST invariance of the gauge invariant operator J and of the vacuum. In the Gribov-Zwanziger approach, 
the situation gets more complicated due to the breaking of the BRST symmetrj@. In the physical limit, £ is no longer a 
BRST invariant operator. In addition, the BRST symmetry is softly broken. Therefore, when turning to physical states, £ 
will no longer be irrelevant, and explicitly influence the value of the correlator. This is not the only observation we can 
make. (x) is not the only renormalization group invariant of dimension 4. Indeed, also the operator £ (x) does not run 
with the scale, as we directly infer from equations ( 1149t and ( 11501 ). We can therefore imagine to study correlators of linear 
combinations of the operators J and £ , where the linear combination is chosen in such a way that the emerging pole 
structure would be real. We notice that this is not a trivial issue in the Gribov-Zwanziger framework Il24l . basically due to 
the fact that the poles of the gluon propagator itself are already not necessarily real-valued. When the Gribov parameter y 2 
is formally set back to zero, we shall recover the correlators of the usual kind in Yang-Mills gauge theories, as the BRST 
symmetry gets restored, as well as the BRST exactness of the operator £ . 

A research project along the previous lines would thus be very interesting to pursue. It would also enable us to show 
that the soft BRST breaking, deeply related to the presence of the Gribov horizon, is not necessarily a negative feature of 
the theory. Rather, it could be very helpful in the construction of suitable operators [39]. We therefore conclude that the 
results in this paper have to be seen as a first step towards the construction of (hopefully) physical correlators in the GZ 
theory. As it should have become clear from this paper, an important tool has been the possibility of embedding the (R)GZ 
theory into the extended model. The nilpotent exact BRST symmetry of the latter model can be used to identify the renor- 
malizable operators by using cohomological techniques, which then also give the renormalizable operators in the physical 
limit. These latter operators will contain the classically gauge invariant operators. At the same time, also renormalizable 
BRST exact operators can be found, which reduce to renormalizable operators in the physical limit, being not necessarily 
BRST exact. It then remains to be seen whether suitable linear combinations of these two types of operators can be found 
that successfully describe physical correlators. This will be the topic of future work. As there are multiple mass scales 
present in the (Refined) Gribov-Zwanziger framework, we expect all of them to influence the pole of the correlators under 
study l39l . 
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